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Abstract 

At the free- fermion point, the six- vertex model with domain wall 
boundary conditions (DWBC) can be related to the Aztec diamond, 
a domino tiling problem. We study the mapping on the level of com- 
plete statistics for general domains and boundary conditions. This is 
obtained by associating to both models a set of non-intersecting lines 
in the Lindstrom-Gessel-Viennot (LGV) scheme. One of the conse- 
quence for DWBC is that the boundaries of the ordered phases are 
described by the Airy process in the thermodynamic limit. 

1 Introduction 

Statistical mechanics models with short range interactions have a free energy 
per unit volume which is independent of boundary conditions. The six- vertex 
model seems to be an exception, at first glance. If one imposes domain wall 
boundary conditions (DWBC), the free energy differs from the one with 
periodic boundary conditions even in the infinite volume limit [7j. 

To explain more precisely we recall the six allowed vertices and their 
weight, see Figure ^ For better visualization only NE and SE oriented lines 
are drawn. The reason for the rotation by tt/4, as compared to the usual 
arrows on Z 2 , will become apparent later. We want to think of the vertices in 
Figure ^as tiles and the equilibrium statistical mechanics as a tiling problem. 
Let us then consider a domain 5* and prescribe the tiles along its boundary 
dS, for an example see Figure E[a). An admissible tiling is such that the 
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Figure 1: The basic vertices of the six- vertex model (first row). By drawing 
only the NE and SE arrows one obtains the tiles in the second row. 



lines are never broken and start, resp. end, only at the boundary dS. In this 
language, DWBC correspond to a square-shaped domain such that the lines 
start regularly spaced at the SW edge and end regularly spaced at the SE 
edge, as in Figure a). The corresponding free energy is computed in 
for finite size corrections see E] • 

In the context of DWBC it has been noted that the six- vertex model maps 
to another, very widely studied tiling problem, namely the Aztec diamond, 
see Figure [2jb). There one has given a diamond shaped domain in Z 2 of size 
2N x 2N, which has to be tiled with 2x1, resp. 1x2, dominos. If every 
allowed tiling has equal weight, for large N there is a central disordered 
zone of linear size N, which is bordered by four regular zones also of linear 
size N and located at the four corners [3]. Thus the limit state for iV — > 
00 "phase" segregates. In contrast, for periodic boundary conditions the 
system is homogeneously disordered, which physically is the reason behind 
the boundary dependent free energies. 

The Aztec diamond has been analyzed in great detail by means of de- 
terminantal processes. In particular one knows the precise statistics of the 
line dividing the ordered from the disordered zone [Hj. The determinan- 
tal property suggests that equal weight in the domino tiling corresponds 
to the DWBC six-vertex model at its free fermion point. The aim of our 
contribution is to study the correspondence between domino tilings and the 
six-vertex model for general domains and boundary conditions, of course in- 
cluding DWBC, on the level of the complete statistics. With this information 
one readily translates the known statistical properties of the Aztec diamond 
to the DWBC six-vertex model at its free fermion point. 

To give an outline: we introduce the domino tiling and the general six- 
vertex model for arbitrary domains and boundary conditions. We explain 
the Lindstrom-Gessel-Viennot (LGV) scheme for nonintersecting paths on 
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Figure 2: (a) A configuration of the 6-vertex model with DWBC The lines 
start at the white dots and end at the black dots, (b) A configuration of the 
Aztec diamond. 



directed graphs and show in which sense the six-vertex model at the free 
fermion point is a special case of the general LGV set-up. Thereby we achieve 
directly the one-to-many mapping from six-vertex tilings to domino tilings. 

2 Six- vertex model with fixed boundary con- 
ditions 

In the usual set-up, at every vertex there are two incoming and two outgoing 
arrows (the so-called ice rule). Thereby one obtains six distinct vertices with 
Boltzmann weights loi, . . . , cue, see Figure Equivalently, one may think of a 
tiling problem with six distinct tiles by drawing only the North-East and the 
South-East arrows. A six-vertex configuration consists in a tiling such that 
the arrows of the vertices match amongst themselves and with the arrows 
prescribed at the boundary of the domain. In terms of the lines of Figure ^ 
each configuration of the six-vertex model can be represented by a set of 
continuous lines which necessarily end at the boundary of the domain. 

In this paper we mainly consider the case where the boundary is a rect- 
angle. Later, we also discuss more general boundary conditions. So, let the 
border of the region we want to tile be a rectangle P$ composed by segments 
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(±1, ±1) between sites of Z 2 , see Figure a). The tiled region S is the one 
enclosed by the rectangle Pg- Define the set dS — Ps H (Z + 1/2) 2 , i.e., the 
mid-points of the basic segments forming the rectangle Ps- Fixed boundary 
conditions means to define the points of dS which are start-, resp. end-points 
of a line. In Figure El^a) the start points are the white dots, the end points 
the black ones. These points are on sides which can be oriented in four dif- 
ferent directions as defined by the direction of the outer normal vector. The 
white dots lie on the SW or NW sides of S, the black ones on the SE or NE 
sides. 

The six- vertex model with fixed boundary conditions has two free param- 
eters. Every elementary tile contains up to 2 line segments of length 1 / \/2 in 
either the NE (/) or SE (\) direction. Denote by iVj the number of vertices 
of type i (the ones with weight Uj) in a given configuration. There are four 
equations reducing the internal free parameters (the weights) from 6 to 2. 
Namely, there are some constants Ci, C 2 , C3, C4 such that 

(1) N 1 + N 2 + N 3 + N i + N 5 + N 6 = Ci, 

(2) 2N 2 + N 5 + N 6 + 2N 3 = number of / = C 2 , 

(3) 2N 2 + N B + N 6 + 2N A = number of \ = C 3 , 

(4) N 5 — Nq = number of /\ — number of \/ = C4. 

The first equation says that the total number of vertices is fixed by S. The 
second and third equation are a direct consequence of the fact that the start- 
and end-points are fixed. Finally, the fourth equation holds because given 
a possible set of lines, one can find any other set of lines by a subsequent 
shift of the corners /\ and \/. During this procedure it is not possible to 
increase (or decrease) the number of \/ by one without forcing the same 
change to the number of /\. 

The weight of a configuration is 

i=l 

and by using the above relations it can be rewritten as 

CupQ* (2.2) 

with C = ^p-( C3+C4 )/ 2 )4 C2+C4 )/ 2 4 C3+C4 )/ 2 ^, u> 2 = uvJt/uwt, and 
0)5 = lo^loq/u^u^. Therefore, w.l.o.g., we can set 

LO\ = Cl>3 = CJ4 = UJq = 1 (2-3) 
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and keep u>2, ^5 as free internal parameters. This choice turns out to be 
convenient for the mapping. 

To a tiling one can associate a height function on S as follows. Fix a 
reference height equal to zero at some point on P$ \ OS. The height along 
the loop Ps is increased (resp. decreased) by one unit when a white or black 
dot on dS is crossed while going to NE or NW direction (resp. SE or SW). 
Inside S the lines of the tiling are the boundaries of domains with the same 
height. In computing the surface free-energy, there are four free parameters 
since two are internal parameters and two correspond to chemical potentials 
determining the average slope, in accordance with our previous argument. 

In the parameter space of the six- vertex model there is a codimension one 
manifold which is determined by the equation 

u>iu 2 + uj-zUJi = uj 5 to e (2.4) 

and is called free-fermion point. Therefore, in case of fixed boundary condi- 
tions, there is only a single internal free parameter. Let us set 

uj 2 = a, lu 5 = 1 + a. (2.5) 

The terminology "free-fermion" refers to a transfer matrix which is the 
exponential of a quadratic fermion hamiltonian. This roughly corresponds 
to line ensembles on a graph conditioned to be non-intersecting (no points in 
common). If we set a = 0, thus 002 = 0, then the free-fermion condition ()2.4j) 
is automatically satisfied, and the corresponding line ensemble fits straight- 
forwardly into the so-called Lindstrom-Gessel-Viennot (LGV) scheme to be 
explained below. For the case a > this is no longer the case and Brak 
and Owczarek [3J approach the problem by introducing "osculating paths" . 
Actually, this is not needed and, as we will explain, the general free-fermion 
case can still be represented in a simple LGV scheme. 

3 Aztec diamond and domino tiling 

A graph closely related to the one presented here appears also in the domino 
tiling of a checkerboard region. A particular example is the Aztec dia- 
mond jSJ. The Aztec diamond consists in a diamond-shaped region D, see 
Figure E(b). The black sites are on even-even and odd-odd sites, the remain- 
ing sites are white. As for the six-vertex model, we can associate a simple 
polygonal loop Pjj enclosing all the lxl squares centered on the elements 
of D. It is a loop which joins the sites of (Z + 1/2) 2 by segments of length 
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Figure 3: Line configurations for the tiling (above) and the line configurations 
(below). The weight of S is a, the weights of N, E, W are 1. 

1. We define the boundary points of D by dD = P D n (Z + 1/2) 2 . For the 
Aztec diamond of parameter N, the set D = is given by 

D N = D NiW U D N>b , (3.1) 

where D^, w and £)jv,& are the white and black sites of the Aztec diamond, 

Dnj> — {(x — m + n, y — n — m) , m — 0, . . . , N , n — 0, . . . , N — 1}, (3.2) 
Dn,w = {(x = m + n, y = n — m — 1), m = 0, . . . ,N — 1, n = 0, . . . , A^}. 

A tiling of the Aztec diamond (the region D^) consists in taking iV(iV+l) 
dominos (2 x 1 or 1 x 2 tiles) to cover every site of Dn (without overlap- 
ping). Depending on the direction of the dominos, vertical or horizonal, and 
their relative positions with respect to the black and white sites, one can 
distinguish 4 types of elementary tiles, called North, South, East, and West 
dominos, see Figure El 

To each of the 4 types of domino one assigns a Boltzmann weight. The 
Aztec diamond has one free parameter. In fact, let N a be the number of 
dominos of type a e {iV, S, E, W} in a configuration of the Aztec diamond. 
Then there are 3 equations reducing the free parameters from 4 to 1, namely 

N(N + 1) = N N + N s + N E + N Wj 

N N = N s , (3.3) 
N E = N w . 

Therefore the Aztec diamond has only one free-parameter. For easy compar- 
ison with the 6-vertex model, we set the weights of the dominos as follows 

wn = we = ww = 1, ws = a. (3.4) 



R 
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Figure 4: The LGV elementary cell of the free-fermion six-vertex model. 



To each configuration of the Aztec diamond one can associate a non- 
intersecting set of lines. This is obtained by replacing the 4 types of dominos 
in the first row of Figure |3] with the lines as in the second row of the same 
illustration. As we will see in the next section, such line ensembles can also 
be associated to the 6- vertex model at its free-fermion point. 

4 Lindstrom-Gessel-Viennot scheme 

Let us first recall the Lindstrom-Gessel-Viennot scheme, see jH] for details, 
and then apply it to the free-fermion six-vertex model. 

Gessel and Viennot start with a graph (V, E) consisting of vertices V 
and directed edges E. The graph has no loops. A path P is a sequence of 
consecutive vertices joined by directed edges. V(u, v) is the set of all paths 
starting at u G V and ending at v G V. The paths P and P' intersect, if 
they have a common vertex. Every edge carries a weight w(e). The weight 
of a path P is 

w(P) = Y[ w(e) (4.1) 
eePn£: 

and we set 

W(u,v)= W ( P )- ( 4 - 2 ) 

P£V(u,v) 

We now consider an r-tuple u = {u±, . . . ,u r } of start-points and an r-tuple 
v = {vi, . . . , v r } of end points. Let Vo(u, v) be the set of all non-intersecting 
r-tuple of paths from u to v. u and v have to be compatible, which means that 
any r-tuple of paths in Vq(u, v) necessarily connects Uj to Vj for j = 1, . . . , r. 
Then the weight of Vq (u, v) is given by 

w(P (u,v)) = det (Wfavj))^.^. (4.3) 
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Figure 5: Mapping from the 6- vertex configurations to LGV line configura- 
tions. The first row contains the 6-vertex elementary tiles, the second row 
the LGV line configurations. Notice that the map from the 6-vertex lines to 
the LGV lines is not one-to-one. We also give the weight of the respective 
elementary tiles. 



For the six-vertex model at its free fermion point, the LGV graph is build 
up from the elementary cells of Figure Thereby, the weights cjj's are as 
required, see also Figure El In fact, the only correspondence which is not one- 
to-one is for the weight u 5 = 1 + a because there are two possible paths, one 
passing by the center of the cell of weight 1 and the second following the edge 
with weight a. Therefore, in terms of non-intersecting line ensembles, we have 
a degeneracy. More precisely, to 6-vertex configurations with iVg vertices of 
type 5, there are 2 Ns different LGV line configurations. However, these are 
only local modifications of each other and any observable not distinguishing 
the precise path inside an elementary cell remains unaffected. 

5 Connection between domino and six- vertex 
models 

The connection between the six-vertex model and domino tilings can be 
seen most directly in terms of line ensembles. Since the mapping from a 
six-vertex configuration to its associate set of lines is one-to-many, also the 
mapping to domino tilings will be one-to-many. In Figure EJa) there is one 
of the 2 5 (since N 5 = 5) possible set of non-intersecting lines associated to 
the six- vertex configuration of Figure EJa). Figure Efb) is the corresponding 
configuration of the Aztec diamond with the same set of non-intersecting 
lines. 
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(a) 



(b) 



Figure 6: (a) One of the possible line-configuration of the 6-vertex model 
with DWBC of Figure Efa). The dotted lines are the edges of the subja- 
cent directed graph of the LGV construction, (b) The corresponding Aztec 
diamond. 



Six-vertex — > domino tilings. 

Given a fixed boundary condition for the six- vertex model on OS, we associate 
the fixed boundary conditions for domino tiling, i.e., we define the domain 
D to be tiled with dominos and the height along the boundary Pp. 

The idea is that the graph LGV supporting the non-intersecting lines 
remains, up to deterministic local modifications, unchanged by the mapping 
to domino tilings. The LGV graph on the black and white sites of the domino 
tiling are 




The region tiled by the six- vertex model and the one of the domino tiling 
cannot be identical. We have to decide whether the squares centered at 
the boundary dS are to be tiled or not. This is achieved by the mapping 
represented in Figure [7| The mapping for the cases E^, k — 1, . . . , 4, follows 
from the fact that, according to the lines for domino tiling introduced in 
Figure El the deleted squares would be part of dominos with a second square 
outside the tiling region. On the other hand, the squares for the cases Ik, 
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Figure 7: The 8 possible boundary configurations of a cell of the six vertex 
model (first row) and its corresponding boundary configuration for domino 
tiling. The dashed lines represents the boundaries. 



k = 1, ... ,4, are part of the tiling region because the second square of the 
domino is already in the bulk of the tiling region. 

The start and end points of the line ensembles are moved accordingly. The 
LGV graph constructed in this way differs from the 6- vertex one by extra half- 
diagonal pieces added or removed (Ei,E^) which have unit weight. 
The other differences between the LGV graphs consists only in edges which, 
due to the boundary conditions, cannot be visited by the lines (Ji, J 4 , E 2 , E 3 ). 
Therefore the allowed line configurations are in one-to-one correspondence 
and have the same weight. 

Particular case: DWBC and Aztec diamond 

A particular case is the six-vertex model in a square domain with domain 
wall boundary conditions, which means that a line starts from each SE point 
of dS and a line ends at each SW point of dS. By the above mapping, the 
resulting domain to be tiled by dominos is the Aztec diamond, see Figure|)]for 
an illustration. This mapping between line ensembles can be easily inverted. 

On more general boundary conditions 

Above we only established the mapping from the six-vertex model in a rect- 
angular region S to a domino tiling. However the construction holds for 
general domains provided the following four boundary configurations do not 
appear. 
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The reason is that the above configurations are mapped into 

and the natural boundary of the tiling no longer coincides with the dashed 
line. 

6 Conclusions and Outlook 

The top LGV line of the Aztec diamond borders, up to errors of size 1, the 
North polar region, i.e., the region containing only North-type dominos which 
is connected with the North border of the Aztec diamond. For size N, let 
us denote the top LGV line by h N (t), t G [0, 2N] (with the axis origin as in 
Figure |2fb)). It follows immediately from the results by Johansson [S| that 
/ijv(^) properly rescaled has limit as a random function. More precisely, let 
a = 1, then 

A n (t) = 2 5 / 6 iV- 1 / 3 (h N (N + riV 2 ^- 1 / 6 ) - N/y/2) + r 2 . (6.1) 
Then An converges weakly to the process Airy process, r > A(t), 

lim A n (t) = A(t). (6.2) 

The same limit process holds for all a > and at all parts of the top LGV 
line dividing the North polar region and the disordered one. Of course, the 
scaling coefficients will change, but the scaling exponent and scaling function 
remain the same. 

The result (16.1)1 . ()6.2|) transcribes to the six- vertex model with domain 
wall boundary conditions. In the geometry of Figure |H1 we may consider the 
outermost six- vertex line denoted by ^dw.tvW- ft borders the region tiled 
errorless with tiles of type 1. Clearly, for all t and N, one has 

\h N (t) -h BWiN (t)\ < 1. (6.3) 

Therefore the limit (|6.1)1 . ()6.2)1 holds also for the free-fermion six- vertex 
model. 
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Figure 8: The phase diagram for the six-vertex model. F are the ferroelectric, 
AF the anti-ferroelectric, and D the disordered phases. The dashed line 
inside the disordered phase is the free fermion line, A = 0. 



The reader might be curious to know about the properties of the outer- 
most line once the free fermion condition is no longer satisfied. For easier 
comparison with the literature, we switch to the standard parametrization 

U>1 = Co>2 = b, Co>3 = UJ4 = O, = UJq = c (6-4) 

with a, b > and exploiting the normalization of the weight we can set c = 1. 
If Figure El we display the bulk phase diagram for periodic boundary condi- 
tions and zero electric fields [HIE]. For DWBC one obtains the same phase 
diagram, however with a very different interpretation. The phase diagram is 
governed by the parameter 

a 2 + b 2 - 1 

A = + °- . (6.5) 
lab v ; 

A = is the free fermion line discussed before. It lies inside the disordered 
phase characterized by — 1 < A < 1. A>lis the ferroelectric phase, while 
A < — 1 is the anti-ferroelectric phase. To understand the structure of typical 
configurations one has to resort to Monte Carlo simulations ^01 Uj • For some 
special points in the phase diagram also analytical results are available. 

In the disordered phase the structure is similar to the one at the free 
fermion line. There are then the perfectly ordered zones at the four corners 
with tiles u>i, U4, u 2 , u> 3 counted clockwise starting from the top and a random 
zone which touches the border of the quadrant at four points. Presumably 
the macroscopic height profile satisfies the Pokrovsky-Talapov law, which 
is a strong indication that the outermost line is still governed by the Airy 
process (see jS] for a discussion). As one crosses to the anti- ferroelectric 
phase, A < — 1, the structure of the outer configurations does not change 
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qualitatively. However, the macroscopic height profile develops a facet at the 
center which consists of a regular arrangements of tiles uj^ and lvq. In contrast 
to the ordered zones at the corner, this facet allows for small statistical errors 
in the tiling. Thus it is less clear whether its (approximate) border ledge is 
still governed by the Airy process. 

In the ferroelectric phase, A > 1, we choose b > a. Then the region below 
the anti-diagonal is filled with tiles of type U2, while the one above with tiles of 
type u>i, separated by a rather sharp interface. The Monte Carlo simulations 
in ^01 indicate that the interface between these two perfectly ordered regions 
has a width which, as N —>■ 00, goes to zero on a macroscopic scale, i.e., the 
width is o(N). If b < a, then the interface is along the diagonal, separating 
tiles of type u 3 from those of type u^. 
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